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Abstract 

Operator product expansion technique is analyzed in abelian and nonabelian field the- 
oretical models with confinement. Special attention is paid to the regimes where nonzero 
virtuality of vacuum fields is felt by external currents. It is stressed that despite the physics 
of confinement is sometimes considered as being caused by "soft" fields, it can exhibit the 
pronounced "hard" effects in OPE. 



1 Introduction 



Inclusion of nonperturbative contributions Q ( proportional to the gauge-invariant local conden- 
sates) in the standard perturbative OPE Q allowed to formulate a powerful method of QCD sum 
rules Q (for reviews see 1^, 0, |5|). Nevertheless some questions in the method were formulated 
by the original authors ^] and still remain unanswered. 

In particular, the relation between the property of confinement and structure of t he sum 
rules series has never been clearly established. On the one hand, one could guess that confinement 
appears as a result of partial summation of some OPE subseries, while, on the other hand, 
confinement itself might introduce some new unconventional terms in OPE series, with the 
structure different from the standard form. 

The phenomenological implication of such new terms, e.g. 0(1/Q^), was investigated in 
1^], where is was related with the short distance nonperturbative physics. The authors of [1lO| , 
11, |l^, checked the role of confinement for QCD sum rules exploiting nonrelativistic solvable 



models, and exact results for Green's functions were compared to the sum rule results. 

Especially popular is the example of nonrelativistic particle in oscillator potential, with the 
Euclidean short-time expansion of Green's function (in 2d, for detailed discussion see cf. the 
3d case in ^) 

Here the first term comes from the free Green's function while the next terms play the role of 
"condensates" namely identifying Borel mass e = one has typical OPE structures: w^/e^ and 

The result (||) has widely been used as an argument that confinement (i.e. long distance 
soft physics) cannot modify the standard OPE and confinement effects should be looked for in 
the partial sums of the type X;~=o Cn(Q^)(-D"F(0)L>"F(0)). 

In what follows we shall demonstrate explicitly that confinement modifies the standard OPE 
for relativistic quark Green's function: new terms appear, which bring unusual power terms in 
OPE. 

It will be shown that the expansion (||) is typical for nonrelativistic potential Green's func- 
tions, while for relativistic particles in the confining fields (or in the confining potential) a specific 
long-distance instability (divergence) occurs in the perturbative expansion, which could lead to 
new power terms. 

Let us stress from the beginning an important difference between OPE in coordinate and 
momentum spaces which was discussed already in original papers jl|, 0] and which will be seen 
clearly in what follows. Studying small x-expansion of a product of two operators (T{ J(0) J(x)}) 
when X — > 0, one observes that in relativistic case (contrary to nonrelativistic one) small value of 
X does not confine virtualities of internal lines in the coresponding diagram in any way. In other 
words, virtual particles created and annihilated by operators J can travel over large distances in 
coordinate space whichever small x is. As a result, in confining theory the product of operators 
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taken at two nearby points carries information about large-distance behaviour of a theory even 
if X is much smaher than typical confinement scale A~^. 

To clarify the mechanism of this phenomenon we start in the next chapter with the Green's 
function of relativistic quark in the linear confining potential of static antiquark, corresponding 
to the Dirac equation with scalar linear potential. We shall expand Green's function in powers of 
string tension (or equivalently in powers of Euclidean time T) and find explicitly a new dominant 
term at small T, and estimate other terms. Comparison with the corresponding nonrelativistic 
Green's function is done and demonstrates that no unusual terms appear in the latter case, the 
expansion being essentially of the same type as in (^). The reason for that is traced to the 
structure of the nonrelativistic free Green's function, for which spacial defiection of particle Ax 
is limited by the time elapsed At, Ax ~ \f^- 

Situation is different however in momentum space. Large external momentum Q plays a 
role of infrared cut-off and if it is much greater than particle mass m and nonperturbative scale 
A, one can successfully perform systematic expansions over m? /Q"^ and A^/Q^. This is how the 
standard OPE technique works. Nevertheless the remaining problem in this case is to determine 
the structure of the latter, nonperturbative subseries. The problem here is that in real QCD 
there are several different nonperturbative scales. The best known are given by nonperturbative 
quark and gluon Q condensates (Vj^) , {F^yF^^y). One can include in analysis higher irreducible 
condensates as well. Another important scale is given by condensate virtualities, see expressions 
(pl|), ( p2D below. So even remaining in the standard OPE framework, one can set oneself the 
task of summation of different subseries in full A^/Q^-expansion. It will be seen below how this 
problem is solved in particular cases. 

Moreover we present a few examples in section 6 where OPE in momentum space starts from 
the terms, which nontrivially account for (monopole) condensate virtuality and hence would be 
considered as subleading in conventional expansion. 

The field-theoretical models are discusses in section 3, where the QCD equations for the 



heavy-light system obtained in the limit of large N(. in |14] are discussed. 



It is shown, in particular, basing on the subsequent results in |15],|16], that exact equations 
have a nonlinear kernel, which at large spacial distances reduces to the linear confining term (j\r\, 
and hence the expansion of the Green's function reduces to the potential example considered in 
section 2. 

We briefly consider abelian models with confinement in section 6 such as QED with monopoles 
and Abelian Higgs model and study influence of conflnement on short-time behaviour of Green's 
functions. We also discuss various approaches related to OPE such as Feynman-Schwinger 
proper time method (section 5) and spectral representations of Green's functions (section 7) 
and study interplay between conflnement and OPE in these frameworks. Finally we present 
short conclusion and outlook. 



3 



2 Relativistic Green's function of a confined quark 



We study in this chapter Green's function of the Dirac equation in the Euchdean space-time 

-i{d + m + a\x\)S{x,y) =6'-^\x-y). (2) 

In what fohows we shah study the function S{x = 0; X4 = 0; y = 0, ?/4 = T) = S{T) as a function 
of T, at smaU values of T. 

The free Green's function S'o(x — y) can be written as 

So{x)= I -^^^Jp:E^{^p + im)=i{m-dm{m^-d^)-^\x) = 



, m Ki{mx) ( X d \ m Ki(mx) , , 

= i{m- d)—^ = i[m t^Htt (3) 

47r2 X V X dx J Att^ x 

where x = \Jx^ + x^ and Ki is the McDonald function. In the massless limit one obtains 

So(-) - (4) 
In the first order in a one obtains in the massless limit 

5(0,r) = So(0,r) +i j d^xSo{Q,x)a\x\SQ{x,T) + ... = So{^,T)+Si{Q,T) + ... (5) 

where function Si can be written in the massless limit as 

ia /■ X I (x — T) 4 



Integration in (pi) yields 



S.(«.r)=p^J^PI^^-Yj4'i^ (15) 



S.(0,T) = £^ (7) 

Consider now the higher-order terms in the expansion (|5|). The typical 0((j"') term looks 

like 

Sn{0,T) = i"" J d'^xi..d'^XnSoa\xi\So...a\xn\So (8) 

It is easy to see that the integrals are infrared divergent at large |x| starting from the term with 
n = 2, however for m ^ this divergence is eliminated and integrals are cut-off by the mass at 
X ~ — . Therefore typical S'„(0,T) has the form for n > 2 



while the n = 1 term obtains corrections of the form 



Sn{0,T)^(^] (9) 



Si{d,T) = ^{mTKi{mT) + 0{mT)). (10) 
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It is instructive to compare (^), (0), (^) with the nonrelativistic expansion (||). One can see that 
apart from difference in free Green's functions, the first dynamical term is nonsingular in the 
nonrelativistic case {§, = while it is singular in relativistic case (0) if T ^ 0. 

To clarify the origin of this difference one can compute 6*1(0, T) for nonrelativistic Green's 
function with linear potential. Note, that free Green's function in 3d is 

a calculation similar to (P) immediately yields 

G-(0,r) = ^ (12) 

which is nonsingular at small T in contrast to 5*1(0, T) in (^). It is easy to see that also all 
higher terms in are nonsingular due to the specific feature of nonrelativistic Green's function 
(pAj): all time intervals are ordered {tn > in-i > in-2) and all space intervals are cut-off by the 
time intervals and the mass, so that quark cannot escape far away during a short time interval 
- in contrast to the relativistic case, when a light quark can travel as far as ^ ^ T for however 
small T. Thus crucial difference between nonrelativistic and relativistic dynamics causes the 
different behaviour of the Green's functions at small distances/times. 



3 Relativistic equation for the heavy-light system 

In this chapter we shall discuss the situation for the field-theoretical model, namely for the 
two-body system made of a spinor particle with the mass m and heavy scalar antiparticle whose 
mass is considered as infinite. We assume that this "meson" interact with confining gauge- field 
background, which is characterized by the Gaussian field strength correlator (see review [41| and 
references therein) 

^£li,M2<^2 = {^r,{F^^^^{zi)<^{zi,Z2)F^^:,.^{z2)<^>{z2,Zi))) = 
lM2'^;^l!^2 ^IJ,iU2^^l2l^l 

) D{zi - Z2) (13) 

where $(x,y) stays for the phase factor 

$(x,y) = Pexp j A^{u)du}^ (14) 

The Green's function of such system can be represented as follows 

{i)[x)^{x, y)i^{y)) = 5o(x, y) + 52(x, y) + .. (15) 
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where 5*0 is given by (|^) while the first nontrivial interaction term has the following form 



TTS2{x,y) 



d'^u J d'^wTi (^So{x,u)iA{u)So{u,w)iA{w)SQ{w,y)^'^ (16) 



where Tr = tr^ tr/, is a product of traces over color and Lorentz indices. We adopt the Fock- 
Schwinger gauge condition with the base point xq = x: A'^{u){u — xq)^ = 0. In this gauge the 
Green's function of the heavy particle is proportional to unity while the gauge field propagator 
takes the form 

{trcA^{u)A^{'w)) = D{0) ■ [{u - x){w - x)5^y - {u - x)y{w - x)^] ■ f{u, w) (17) 
where dimensionless function f{u,w) is given by the following expression 



fiu,w) 



1 



D{0) 



1 1 
daa / dP(3D{au — (3w) 



(18) 











Functions of this kind are often used in the formalism of coordinate gauges, one can find in 
Appendix of the present paper detailed analysis of f{u,w) for particular choice of Gaussian 
ansatz D{z) = L'(O) exp(— z^/T^). We shall keep only the function D{z) in what follows since 
the function Di{z) is not responsible for confinement effects. It was also found on the lattice 
that nonperturbative part of -Di(z) is significantly smaller than that of D(z) in QCD, see 
and references therein. 



In momentum space ( |16D takes the form: 



/d'^l f d'^k f f 
J^:^ J J^:^ J J l^^f^ exp{il{x - y))- 



d d 
drp dsa 



r — ka 
s — kl3 



{I — s + rY + rn^ {I — s)^ + m'^ 



•((5pcr(3m^ - - 2lr - sr + s^) + 4lpl„ - ilpS^ + 2lpra^ - 2l„Sp + 2spSa^ - r^Sp - VpS^ 



(19) 



The properties of the expression (|T^ are determined by the interplay of external parameters 
such as particle mass m and distance |x — y| and properties of the confining background encoded 
in the function D{z). In case of QCD the latter is usually found on the lattice |45]. It decays 
with distance and has some typical correlation length scale which we denote as Tg throughout 
the paper. The exact dependence of D(z) on z is of no principal importance, one usually 
takes exponential fits (see [^]). At the origin D{z) is normalized to the nonperturbative gluon 
condensate, according to 

^(0) = ^ {treF^.F^,) 

It is worth mentioning that the actual numerical value of Tg in gluodynamics and QCD is rather 
small: it is estimated as 0.22 Fm for quenched SU (3) and as 0.34 Fm for full QCD with four 
flavours |45, 43, 41|. As it will be clear from what follows this circumstance bounds region of 
applicability of conventional OPE based on local condensates. 
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We study first the heavy quark case, i.e. we assume that mTg ^ 1. The integrals in (|16|), 
( |T9| ) are saturated at momenta P of the order of the mass m? and one can make systematic ex- 
pansion over 1/mTg. Straightforward although rather lengthy calculation leads to the following 
answer for the heavy quark condensate: 



TV52(a;,a;) 



247r2 



m 



44 1 
45m2f2 



+ 



1 



(20) 



where Tg is defined as 



7-'2 
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4D{0) J (27r4) 



(ticiFD^F)) 
(tr,F2) 



(21) 



where the last relation is valid in Gaussian approximation when all contributions from higher 
correlators are neglected. Let us mention that virtuality of quark condensate usually measured 
by the quantity 



A2 



(22) 



in the sum rule approach is comparable with that of the gluon condensate (|2l|), indeed, = 
0.4±0.1GeV2 according to El, while Tg was found on the lattice to be 0.34±0.02Fm in SU{3) 
with 4 dynamical flavours [^ ], i.e. XqTg is of the order of one. It could be instructive therefore 
to reconcile our approach with the method of nonlocal quark condensates worked out in 



For D(z) oc exp{—z'^/Tg) with the correlation length Tg (as is used in Appendix) one has 
Tg = V^Tg. The first term in the expansion (^) is well known OPE-result for the heavy quark 



condensate 1^^, see also |4^. The second term is the first nonlocal correction. It is worth 
mentioning that due to the smallness of Tg (see above) it can be omitted as numerically small 
correction for b, t quarks only, while for s, c quarks keeping only the first term in the expansion 
in 1/mTg is not to be considered as good approximation. 



Equations of the form (|21|), (|22| ) account for nonzero virtuality of vacuum lines in standard 
OPE language - one considers quantum averages, which contain derivatives. As we shall see in 
what follows, this language is not universal and implicitely assumes small averaged virtuality 
corresponding to the vacuum state, i.e. large Tg limit. Another essential ingredient of this 
language is the use of equations of motion for such averages. Although it is rather easy to 
justify the validity of this component of the approach in abelian case, to the best of author's 
knowledge, this procedure has never been proved for nonabelian theories with the level of rigour 
adopted in the field. Since we are discussing nonlocal correlators, the following remark is of 
importance. Consider parallel transported field strength tensor -F^jy, i.e. 

G^uix,xo) = ^{xo,x)F^^{x)^{x,xo) 

and nonlocal gauge-invariant two-point correlator 



{trcG{x,xo)G{y,xo)) 



(23) 

The above correlator depends on the positions of the points x,y,XQ and on profiles of the 
contours used in factors <I>. However, if x — > y all these dependences disappear (phase factors 
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cancel each other, while x-dependence is prohibited by translational invariance) and the resulting 
local average coincide with {tVcF"^). Let us consider now expansion of ( p3| ) if \x — y\ is small. In 
principle one might consider two different expansions, with correlators involving derivatives in 
both cases. In the first case it reads: 



{ticG^ij{x, xo)Gpaiy, xq)) w (trc-F^) + {y -x)a- { trcG^u{x, xq) 



where the derivative is given by 



dGpa{y,xo) 
dya 



+ 



(24) 



y=x/ 



^%^^ = ^(xo,y) (^D^Fp^{y) + iiy-xo)p J sds[GpUz,y), F^^iy)]^ $(y,xo) 



(25) 



and [..,..] in ( pq ) denotes commutator with respect to the color indices. The second term (and 
all higher terms) in the r.h.s. of (|2^ ) contains nonlocal part and depends on contour profiles and 
on the position of the points x, xq unless x = y = xq. On the other hand the expansion goes in 
powers of the quantity (y — x) which is assumed to be small. 



In the second case one expands each G in (|23| ) in the vicinity of the point xq: 

{x,XQ)Gp^{y,XQ)) 

K. {iTc{F^i,{xQ) + {x - Xo)aDaFpuixo) + ..){Fp„{xo) + {y - Xo)aDaFp^{xo) + ..)) ~ 

« (trcF^) + {y - Xo)a ■ {tTcFf,uixo)DaFpa{xo)) + (x - Xo)a ■ {t'CcDaFp^{xo)Fp^{xo)) + ... (26) 

This is an expansion adopted in conventional OPE. In contrast with (^) nonlocal parts are 
absent, the price to pay however is that the expansion goes in x — xq , y — xq instead of y — x. 
Needless to say that in many physical applications |y — x| can be small whereas \x — xo\ and 
\y — xq\ are very large. Notice also that opposite situation is impossible: smallness of |x — xo|, 
|y — xol implies smallness of \y — x\. 

After this rather academic discussion we come back to the limit of small quark mass and/or 



correlation length mTg <C 1, which is opposite to what has been explored in (|20|). As it was al- 
ready mentioned, in real QCD parameter (trc-F^)T^ can be considered as small, even in presence 



of dynamical quarks. In particular, typical momenta P in ( |l9|) can be rather large in comparison 
with nonperturbative scale given by the condensate y^^tivF^ but still small when compared to 
nonlocality scale 
regime. 



Tg ^. Test particle resolves nonlocality of vacuum field correlations in this 



The Green's function S2{x,y) we are interested in is defined in ([iq). We are working in 
coordinate representation here and choose the Fock-Schwinger gauge reference point xq at the 
origin xq = x = 0. We rewrite S2{x = 0,y) using (^) as 



S2{0,y) 



iPjO) 
647r6 



d u d w 



— -2—]- {4{uw) -uw) 



■f{u,w) 



u — w 



+ 



m 



{u — w)^ (u — w)"^ J \ {w — y)^ {w — yY 



w - y 



+ 



m 



(27) 
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where we have kept only hnear in mass m terms in propagators since we consider small mass 
limit. The kernel f{u,w) is defined in ([T8|) and u = u^j^. 

The actual value of this integral is defined by the properties of the function /(u, w) which 
encodes all nonperturbative dynamics in the chosen Gaussian approximation. They are rather 
peculiar however (see Appendix) and this circumstance precludes one to obtain exact analytic 
answer. On the other hand, ( p7[ ) can be calculated numerically for any particular ansatz for 
D{x). Let us investigate general structure of 82- In massless limit one immediately obtains 
lim 5*2(0,1/) = according to absence of chiral symmetry breaking in the problem in question. 

It is seen that ^2 is UV-finite (small n, w domain) because nonperturbative background is soft: 
lim f(u,w) = 1/4. In infrared domain ^ Tg the integral is convergent due to the 

U,W—rO 

properties of f{u,w) (see Appendix). One obtains in massless case the following leading term 
at small \y\: 

S2{0,y) = -ic-D{0)-y + O{y^) 

Numerical constant c is determined by the function /(u, w), but is Tg-independent. The massive 
parts of 52 provide finite contribution at y = 0: 

S2{0,y) r.imD{0)T^ 

If mass is increasing and reaches values of the order of T~^, it begins to work as IR cutoff instead 
of Tg and one comes back to (|20| ) . However if mass is small then light quarks essentially feel the 
virtuality distribution of vacuum gluon fields (i.e. the profile of f{u,w)). 

It is instructive to show how the potential problem considered in section 2 appears from field- 
theoretical framework invoking by us here. To this end one is to consider equation for heavy-light 
system which was obtained from the QCD Lagrangian in [14] in the limit of large Nc- Keeping 
only the Gaussian field correlator one has instead of (|2|) the equation for the quark Green's 
function (made gauge-invariant due to phase factor coming from the heavy source propagator) 

-i{d + m)S{x, y)-i I M{x, z)d^zS{z, y) = 5^^'^ (x - y) (28) 



where the nonlocal kernel M{x,z) depends on the exact Green's function S{x,z), making Eq. 
( p8| ) nonlinear. Till the end of this section we are working in the so called modified Fock- 
Schwinger gauge (see all details in pH]) where the temporal axis is singled out. We have retained 



for simplicity only color electric part of the correlator as defined in |21| {Ei{x)Ek{z)) ~ 5ikD{x — 
z). Assuming for D{x — z) Gaussian ansatz, one arrives at the following form of nonlocal kernel 

M(x,y), 

M{x,y) = Dmxy)f{^.y)S{x,y)e^v (^-^^^^^^ (29) 

where /(x, y) is given in Appendix. Notice that x, y are three-dimensional vectors here and not 
four-dimensional as in (|2T 



As it was shown in [15| using the relativistic WKB method developed in |14|, the function 
S{x,y) at large x,y, i.e. if \y\ ^> Tg can be written in the following form 

Sih, f , y) = ie-(-\'\+^^'^gix, y) ( ^^^^ \ (30) 
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where h = X4 — and g{x, y) is a smeared J-function 

(7(rc,y) =p)(x-y), af^ ~ ^ry^ > 1, (31) 

moreover for large x and y, and |x — y| <C |x| ~ |y| (see and also Appendix of the present 
paper) 

(32) 

\x\ 



and in this region one can integrate in ( p8| ) over d^z, since 

M(x,z) ^CT|f|(^(^)(a;-z) (33) 



where the string tension a = {tt /2)D{Q)Tg for Gaussian ansatz (see (32) below). Thus at 
large spacial arguments the kernel M coincides with linear potential considered in the previous 
chapter. Therefore all estimates for terms in the expansion proportional to M'^,n > 2 are in 
agreement with those for the local case, eq. (9), since integrals in these terms are essentially 
saturated by large spacial distances, l^^"^! » Tg. 



Although potential behaviour (|3^) is typical for large-T - regime, it is instructive to show 
how the nonlocality cures l/T behaviour found in local potential problem. We shall demonstrate 
now that Green's functions in question have finite limit when T — > either for small or for large 
Tg. Let us briefly analyse the situation with the nonlocal equivalent of (^), i.e. 

s!"'>(o,T) = 11/ d*. / <i* /''''''"j^/frf"''''''"' °''P {-^) <^^' 



with S given in (|3y). It is convenient to introduce dimensionless quantities 

x^i = x^^/a, y^ = y^\/a, fg = Tg^Ja, f = T^/a. (35) 

Rewriting (^) in terms of tilde variables, one immediately realizes, that S tends to 3-dimensional 
delta-function only in the limit when |y| ^ 1 and otherwise the integral is defined by the 
region \x\ ~ \y\ ~ 1, when the nonlocality is at work, i.e. |x — y| ~ \y\. Imposing the limit of 
small Tg, i.e. <C 1 one reduces two powers of x,y in the numerator of (|3^), but the integral 
is still defined by large values of x,y of the order of unity and one finally obtains 

5f^^(0,r) -const -(73/2 (36) 

One can also show that the same estimate holds true also for higher terms 0(m") and write 
Sn^^\o,T) ~ c„cj^/2. Consider now the opposite limit fg » 1, i.e. Tg » 1/v^- this case fg 
does not confine the differences x — y in f{x,y) and X4 — y^ in the exponent in ( p^ to small 
values as compared to |5;4| or |y|, |y4|. Therefore the integration over d{xi — y^) is limited 
only by the exponent in (|30|). As a result one obtais for Si the following estimate (as always, 
we assume mass ra to be not large, m < ^/a) 



5;^)(0,r)~ const (37) 

Thus in both cases the normal procedure of OPE, based on the analysis of subsequent terms of 
perturbative expansion (with separating soft and hard parts of diagrams) is not applicable and 
one should sum up the whole series or else solve the nonlinear equation (^) exactly. 
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4 Another field-theoretical example: how the linear confine- 
ment is built up out of condensates 



In this chapter we consider another example: a scalar (Higgs-type) particle interaction with the 
background Yang-Mills field, and again calculate the heavy-light Green's function, where light 
particle is the color fundamental Higgs, while heavy source is like an antiquark. The Lagrangian 
and the Green's function are given by 

L = ]{Fp' + \D,^\' - G(^)(0,r) = immTMT)) (38) 



One can rewrite G^'^^ as 

G(^)(0,r) = ((m^ - Dl),},mT))B (39) 

where we have introduced the external (background) field B^: = — igB^. For simplicity 
of consideration we shall confine ourselves to the vertices 

U^g\Bf^'){B-^^T (40) 
and choose the gauge to write equation for G^'^\x,y) analogous to (p^), 

(m^ - dl)G^^^ (x, y)+ f I{x, z)G^^^ {z, y)(fz = 5^^) {x - y) (41) 



where I{x,z) I^'^\x, z), and I^^^ refers to the kernel with one power of B^, while 

/(2) corresponds to the Lagrangian ( ^0[ ) and can be written as 

I^^\x,y)=6^'\x-y)g'Blix) (42) 

The contributions from I^^^ not analyzed by us here are of the same general structure as that of 
I^"^^ apart from nonlocality controlled by the particle mass m. Our consideration in this section 
is of illustrative purpose only. 

Let us take now the first order graph in g"^ (keeping only /^^^ in (|4l|)). In Euclidean space- 
time one has 

G'f\0,T)=g'Q d^xG^^\0,x)Bl{x)Gi^\x,T)'j ^ (43) 



where 

47r2 



g''^\x) = -^-Ki{mx), X = Vx^. (44) 



We shall be interested in the vacuum averaged expression for G^^ and to this end one 
should express B^ (x) in terms of field correlators (one way) or in terms of condensates (another 
— standard way). In the gauge [^] one can write e.g. for B|: 



X 



{Bi{x))= du^ dv,{F^4iuMu,v)F,4{v)) (45) 
Jo Jo 

and using [21| and keeping as in section 3 only the confining part D, one has 

(Blix)) = dui dviD{u-v), i = 1,2,3. (46) 
Jo Jo 
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By way of example let us consider exponential ansatz for D[u — v). As it was already said, 
this behaviour of D{x) was observed in lattice simulations at distances larger than some typical 
correlation length Tg. So one has 



{bI{x)) = D{Q)x^ da d(3ex];) 



\a-f3\ 



(47) 



Notice absence of additional multipliers a , /? in ( p7[ ) contrary to (^), it is a property of modified 
Fock-Schwinger gauge ||2^ (temporal axis is singled out) used by us in this section, instead of 
usual one, used in ([l^) (one point is singled out). 



(48) 



Straightforward integration leads to the following result: 



{Bi{x))=2D{0)\x\Tg 1 



11. 

\x\ 



1 — exp 



At large \x\ ^ Tg, one has from 



{Bjix))^2D{0)Tgi\x\-Tg 



(49) 



Notice that nonlocality enters ( p9[ ) in explicit way. At small when \x\/Tg ~ 1 the linear 
behaviour of (|49|) is replaced by the quadratic one 



(Biix)) ^ D{0) ■ X 



^2 



(50) 



It is clear that vacuum field correlator method implemented in (|45|)-([49|) demonstrates the cre- 
ation of the string between the Higgs particle at the point (x, X4) and static source at the point 

(0,X4). 

Now let us look at the same problem from the point of view of standard OPE. According 
to general rules [Q], one should expand F^j^{u),Fy^{v) in (|45| ) in the vicinity of a point 

-u = i; = 0, U4 = t;4 = X4 (in the Fock-Schwinger gauge that would be the point usually chosen 
at the origin, zq = 0). In this way one obtains 



kBi{x)) = y: 



1 



{n + iy.{m + iy. 



XiXi^...Xi^XkXk^...Xk^{Di^...Di^Fi4^{0)Dk^...Dk^Fki{0)) (51) 



In this form ( [5l| ) the appearance of the string is not visible, and one should rearrange the 
derivatives in nontrivial way, so as to separate out the correlator D{u — v) as in (^). Derivatives 
of the latter produce powers of T~^, while dependence on the sum i(nj + Vi) in the integral ( |4^ ) 
is separated out to yield linear confinement in (|49|). 

Now we consider the expansion of G(^) in powers of . From (^) and (|4^), (Q) it is clear 
that one obtains 

gS-)(0,T)~5^^, G^^^ - {^X m\ (52) 
m \m-^ J 

where nonlocal constant C4 ~ D{0)T3 ~ f dzD{z). All integrals like ( |4^ ) are diverging at large 
distances for m = and cut-off at x ~ when m 7^ 0. 
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It is instructive to turn to the momentum space and define the fohowing one-dimensional 
Fourier-transformed Green's function: 



G^^\q)= J dTG'f'\0,T)eM-iQT) (53) 

— oo 

(Our problem is 3-1-1 dimensional, we do not perform 4-dimensional transformation however 
since the temporal axis was separated by our gauge choice from the beginning). Since {B1{x)) 
does not depend on the temporal coordinate x^, integration in (Esl) is trivial: 



<'(M. . ^/.^.^^E™ (.|,„„ ^ ^ ^ . ^ (54, 
where = Q"^ + m?- Expression ( |54|) has the following asymptotic expansions: 

o2 



and 



GfiM) = I- • (1 - 2MT3 + 4M2r2 + 0{M^T^^) , MTg < I (56) 

It is clearly seen that the actual answer is given by different series in regions QTg <^ 1 and 
QTg » 1 (we assume that Q ^ m and M ^ Q). The expansion (^5|) is associated with 
the standard OPE (]5l|), while (^6|) goes essentially in powers of nonlocal quantity. It is also 



worth mentioning that both expansions ( p5|) and (56) are model dependent beyong the leading 



condensate term and actual coefficients in (pSj), (pq) are determined by the profile of D{z). 



5 Feynman-Schwinger formalism and OPE 

We start with the same Green's function G^*^^ (0, T) and write the Feynman-Schwinger represen- 
tation (FSR) for it (see |22],|23] where refs. to earlier papers are given, for more discussion see 
12 



G(^)(0,r) = dsj {Vz)o,TeM-K){W{C,)). (57) 

Here K = m?s + 5/0 ^"^dr, s is Schwinger proper time and {W{Cz)) is a Wilson loop consisting 
of a straight line (0, T) and the trajectory of the Higgs particle from to T. Notice also, that 

(^^)o. = n ^ / 1^ »P ('P (E ^4n) - t) ) (58) 

In Gaussian approximation {W{Cz)) can be written as 

{W{Cz)) = exp (^-^ J^da^,{u) J^dapx{v){F^u{u)Fpx{v))^ (59) 

and we have omitted for simplicity the parallel transporters inside {FF). Here S is the pre- 
scribed minimal area surface in the loop Cz (there is no sensitivity on the choice of S when 
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all higher cumulants are kept; with the choice of the Gaussian correlators and minimal surface, 
the contribution of all higher correlators was estimated to be around few percents, see |41] and 
references therein). 

For small contour Cz (which means that not only T is small but also spatial size of the 
contour is small), one has from ( |59|) 



{w{c,)) = exp ( -|!^(f;,(o)f;,(o)) 1 (eo) 



For a rectangular contour Cz of an arbitrary size R x T one can write 

{W{C)) = exp (^-^J d^x <fyD{x - y)^ (61) 
Choosing for simplicity D{z) = L'(O) exp(— z^/T^)), one has 



1 /'n..^^2._v^D(0)^2 



2 



D{z)d'z = —T^ (62) 



and finally 



where we have defined 



with the expansions 



(iV(C))=exp(-^Lg)L(S]] (63) 



oo 



L{u)= I die-*'"^ (64) 



Liu) = xf^{l-T^ + ^ + ...}, u»l (65) 



72^2 80^4 



L{u) =TT + 0{^), u<l. (66) 

Now we consider (W) inside the integral (^). If one assumes, that for small T one can 
indeed use the approximation of small area of the loop Cz, i.e. Eq. (|60|), then one has in the 
relativistic case, but considering T small, T <C — and expanding 



C^^\0,T)=C^^\0,T) + C^;^\0,T) + ... (67) 

where 

g[^\0,T) ~ g\S^){F''F'')^ ~ g^T^iF^F''). (68) 

(we have assumed according to what was said before (|60| ) that (5^) ~ T^. This is a standard 
result of OPE analysis with a local condensate accompanied by higher powers of T (or higher 
powers of in the momentum representation). 

Let us now consider again the term 0{g'^), but now taking into account the nonlocal char- 
acter of the correlators {F{x)F{y)). To this end we expand the Wilson loop in ( |57| ) and making 
use of two simple identities 
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/•oo rs rT\ poo roo roo 

/ ds dn dT2f{s,Ti,T2)= ds dTi dT2f{s + Ti+T2,Ti+T2,T2) (69) 

JO Jo Jo Jo Jo Jo 

one can write 



Gr{0,T) = j d'^xj d^yGr(0,x)(^^(x)yl,(y))x^G^^^(x,y)y,G[,^^(y,T) (70) 

Notation used in (|70|) implies, that x^{t) = and {A^{x)Ay{y)) is expressed through a 
vacuum average of field strength {F{u)F{v)), e.g. as in the coordinate gauge 



{Ai{x)A4{y)) = r du, r dvk{E,i{u)Fk4{v)). (71) 
Jo Jo 

One can show (see, for example, Appendix B of [^) that /dx^ and one recoveres in 

(|70| ) the usual perturbation expansion for G'-'^'^ , where now in contrast to the chapter 4 only the 

linear vertices '^'^A^ d ^ are taken into account (the term (f'^A'j^ would also appear in (57) when 
one takes into account term with x = y). 



From ( |70| ) one can deduce that the r.h.s. stays constant at large |x — y|, while it decreasing 
for large \x4 — 7/4I (this is especially clear when one uses for the correlator D(u — v) the Gaussian 
form). Therefore the integral (^) is convergent both at large x, y, and at small x, y. Integrating 
( |70| ) one obtains 

G[^\0,T)r^T^{gF'^gF'^)r^a, (72) 

since -D(O) ~ {gF'^gF"') ~ '^/^g- Thus one obtains a nonlocal constant for small T <^ Tg. 
Comparing ( |68| ) and ( [72D one can see that at small T the correct (nonlocal) procedure yields 
a larger (dominant) term as compared with the standard OPE estimation. The reason again 
lies in the fact that relativistic trajectories occupy larger area for the Wilson loop when treated 
perturbatively and nonlocally, whereas in standard OPE treatment one attributes to this term 
the local condensate, implying that the Higgs particle (or quark) does not go far from the static 
source. 



6 Remarks on OPE in abelian theories with confinement 



We are going to discuss OPE in abelian confining models in this section. The complications 
due to path ordering are absent in abelian case and one may consider general expression for 



the two-point correlator of the field strengths in the form (13) where the functions D{z), Di{z) 
depend entirely on z = x — y. We assume that confining properties of the theory are caused by 
condensate of monopoles, hence the equations of motion take the form: 

d^F^iy = ju ; df^F^y = (73) 

where F^y = \e^j,ypaFp^ and j^, are electrically and magnetically charged currents, respec- 
tively. We define polarization operator n(g^) of the electric currents as 

d^x{jp{Qi)j^{x)) exp{iqx) = (^9^ - Qt^quMq^) (74) 
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Differentiating abelian analog of (13) and taking into account equations of motion, it is straight- 



forward to obtain the fohowing relation: 

n(g2) = J d^x (d{x) + ^Z)i(x) + 2;'^) exp(iQx) (75) 

In d = 4 case it can be rewritten in symmetric form as 

{j^U^)) + {JMM^)) = -\ - 5^a.)(F,,s(0)F„^(x)) (76) 

where for the condensate one has 

{Fo.p{^)F^p{x)) = 6{D{x) + D{x)) (77) 

The function D(x) = D{x) + 2Di {x) + x^^^ corresponds to the confining part of the correlator 
of dual field strengths F in the same way as D{x) corresponds to the correlator of F in (0). In 
case of d = 4 QED without monopoles, one easily finds 



, , 1 /e2(x) 1 de2(x2)\ 



For polarization operator one obtains 



nC'?') = / d'^^Mn^)^^^ (79) 



It is evident that free field term ~ Cg/x^ does not contribute to n(g^) and the only nonzero 
contributions to the r.h.s. of (^) comes from either the running of the charge e^(x^) (as in 
perturbation theory) or from nonperturbative parts of D{x),Di{x)^ if they are nonzero. Let 
us examine the latter contributions to H{q^). Standard OPE reasoning would suggest to look 
for leading term of this kind in the form {F'^)/q^. It is easy to see that for functions D{x), 
which are smooth at the origin (for example, D{x) = -D(O) exp(— x^/T^)), the corresponding 
contribution to n(g^) is exponentially suppressed at large (i.e. for q'^Tg ^ 1), it means that 
power corrections are absent in this case, in other words nonperturbative background is "too 
soft". In particular, there is no D{0)/q^ term. For D{x) such that it is not smooth but finite at 
the origin (e.g. for often used exponential fit D{x) = D{0) exp{—\x\/Tg), one has as a leading 
large-g nonperturbative asymptotics 

The situation becomes even more dramatic if D{x) is singular at the origin (as it happens, for 
example, in the London limit of Abelian Higgs model Q]), where D{x) ~ M^/x^ if x ^ 0, one 
has in this particular case 

^n(g2) = oc ^ then g2»M2 (81) 

This 1/g^ regime in AHM is bounded from above, however, by the Higgs mass rriH- if q'^ ^ "t-h? 
the Ginzburg-Landau description of the condensate is not valid, broken symmetry is restored 
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and microscopic degrees of freedom come into play. Presumably the same reasoning in applicable 



to "thin" strings scenario, proposed in |46]: at distances much smaller than coherence length 
neither "thick" nor "thin" strings can be formed. Notice that the string tension a depends on 
mu logarithmically in the London limit: a ~ \og{mH /M). 



It is interesting to compare the result (|8lD with an answer for massive photon propagator. 
It can be obtained from ( |7^ ) taking e^(x) = mxKi{mx) which corresponds to massive vector 
field propagator {Afj_(0)A^(x)) . Differentiation in (|79|) yields 



m 



2 



U(q^) oc TT then > (82) 

This result is obvious from the form of propagator in momentum space. Notice the sign difference 
between ( p2| ) and (|8T[). It can be said, following ||3^ that leading power correction AP[q'^) in 
confining theory is caused by exchange of massive particle with tachyonic mass. This interesting 
point will be discussed elsewhere. 



7 Spectral representation of Green's functions and OPE 



In this chapter we shall look at OPE from another point of view, trying to calculate terms of 
OPE using the known properties of spectrum of gauge-invariant Green's function. 

This type of analysis was done most extensively for the 'tHooft model (1+1 QCD at large 
Nc) where exact results for the spectra and Green's functions are known. (For details of 



analysis the reader is referred to 28, 29, 30, We follow most closely notation and the 
line of reasoning of |31]. We consider again the heavy- light system but now in the d = 1 + 1. 
The Green's function can be written as 



oo 1 

G(««)(x) = (0|g-(x)$(x,0)g(0)|0) = Y: 7^(fe^M)'^> 



n=0 



Defining on the other hand the correlation function 

P{q^) = i [ e"^^d'^x{0\T{qQ{x),Qq{0)}\0) i I G^^^\x)e''^''d^x, 



one can write and expansion in inverse powers oi E = niQ — qo 



+ pert. part 



(83) 



(84) 



(85) 



where Pq = iDq. 

On the other hand one can write a spectral decomposition (dispersion relation) for P{E) 



P{E) = ^moV^Y: 



1 



E + En vr '^^/n{,/n + e) 



(86) 
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where we have used notation mg = g^N^-jix^ 2mQ^pTie = E and relations for the heavy-light 
spectrum [^] 

= 2m,V^ (l + O i^^)) , fl = V^n{l + {^^)) . (87) 
Now one can compare (^) and ( |86|) and expanding the latter in powers of ~ pr, one 



obtains |31] for coefficients in 

{qP^^q) ~ (g-g)(7rmo)2"n! 

The factorial growth of coefficients in ( p8| ) is typical both for 1+1 and 3+1, as will be shown 
below in this chapter. 

One can do another derivation of the coefficients (|8^ ) starting from equations of motion in 
which case instead of (^) one obtains 

2/ 



{q{x,D,f"q) ~ x-'^nliqq) f"^^ 



Thus appears another feature (or a puzzle, as it was formulated in |31]): condensates 



computed from the spectrum or from microscopic equations of motion have drastically different 

scales: tjiq in the first case and [j^j in the second case, where mg tends to zero in the chiral 
limit. 

We shall now show that in the 3+1 QCD at least for Nc ^ oo the situation is very similar 
to that of the 'tHooft model: 

a) OPE coefficients of the expansion ("condensates") have factorially growing be- 
haviour. 

b) Condensates calculated from spectrum and from diagrams (plus equations of motion) 
are different. 

Consider now the 3+1 problem - description of the selfenergy part Il{q'^). For two light 
quarks the standard OPE of Il{Q^) in the Euclidean region is well known |||] 

^,^2n 1 / as\, Q2 6m2 2m(qq) as(FF) , , 



Following |32] one can use the background perturbation theory for the calculation of n(Q ) 



and represent it in the form 

n(Q2) = nW(Q2) + a,n«(Q2) + alXl^''\Q\ (91) 

Let us first consider n('')((5^) (for details of computations the reader is referred to |3^ and 
papers quoted therein). 
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In the large Nc limit II^^\Q'^) has the form 

n=0 ^ 

The masses M„ can be taken as the eigenvalues of the well-known Hamiltonian, which was 
derived from QCD with the assumption of area law for minimal surface and was shown to be 



valid for small angular momentum L = 0, 1, 2 [33|, while for larger L the string rotation should be 



taken into account, AHstr, yielding the correct Regge slope {2-Ka) ^ for masses M„ 34, R^] 



H^^^^n = MniJn] H^^^ = 2.'^ + m] + ar + AHstr (93) 



Solutions to (^) can be written in the form 

= 27rCT(2n^ + L) + (94) 

where Mq ~ m^. For C„ one has 

Cn{L = 0) = '^Q}N,ml Cn{L = 2) = ^Q}N,ml (95) 



Here ttiq = 47ro", and Qj is the electric charge of quark of flavour /. Taking into account 
degeneracy of masses with L = 0, n,. = 1 and L = 2, n,. = the total Cn is the sum 

Cn = Cn{L = 0) + Cn{L = 2) = Q}N,ml (96) 
Since Cn does not depend on n in this approximation, one obtains the sum 

E 772^ = ( ^-^] + const (97) 

where the constant term is divergent and is eliminated by renormalization of Il{Q'^) — > Il{Q'^) — 

n(o). 

Here il^{z) is the Euler function 

^ ^ fc=i 

where -B„ are Bernoulli numbers. Hence at large the leading term in (^) yields 



For S> Mq this term coincides with the leading term in the OPE (^0[) (the latter is written 
for Qf = 1). 
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From (^) and ( |97| ) one can compute also the next terms of the expansion in ^ 

^ n=l ^ 

It is clear that A„ at large n grow factorially due to the asymptotics of Bernoulli numbers, 
Two properties are clearly visible in the expansion ( |100| ) 

a) the "condensates" ttiq" are large, mo w 2.5 GeV, as compared to the standard OPE 
condensates, e.g. {FF) ~ 0.1 ^ 0.2 GeV^ 

b) the coefficients A„ grow factorially, which is in agreement with discussion in |2£] and 



analysis of the 'tHooft model in |3C, 31|, signifying that the OPE series is asymptotic. 



Thus in both cases, 1+1 and 3+1, when confinement is present and spectrum contains 
nondecreasing probabilities C„ (which is the feature of linear confining interaction) the OPE is 
a factorially diverging series, implying renormalon singularities in the Borel plane [^]. Another 
feature which is general to both 1+1 and 3+1 theories, is the mismatch between condensates 
calculated via spectrum (as in ( [lOOD and via diagrammatic analysis (as in (pO|)). In |31| a 



possible solution of this masmatch for the 1+1 case was suggested, which introduces the notion 
of "effective condensates", which may differ from actual condensates (defined, for example, on 
the lattice) due to the asymptotic character of the OPE series. 

In 3+1 case there is another possibility to explain the mismatch, namely one should take 
into account that coefficients of all higher condensates get contribution not only of the leading 
terms in n of M„ and c^, but also subleading terms, and the final result for say A4 could be 
two orders of magnitude smaller due to cancellation between different terms, thus removing the 
mismatch. However this requires a mechanism of fine tuning between the subleading coefficients, 
the physical reason for which is still not known. 

One could leave discussion of the mismatch at this point, if another check were not possible. 
Indeed, let us take the OPE with large (spectral) condensates and do a sum rule analysis of 
experimental data for e'^e'^ hadrons with 1=1 (see 



This analysis was done in [^2| using the hadronic ratio R^{s) = 127rImn^(S'). For 1 = 1 
adding the perturbative terms with known coefficients as in [||, but taking the background 



modified coupling constant ||3^ e.g. in one loop 

asiQ^) = — 



An 



bo In 



"AT 



where Mb ~ l.SGeV, A^^ « 482MeV, one has 



= I E Ct^Sis -M^J + Ul + ^ + 1.64 ( ^ ) ) (101) 



ra=0 



C^-^=m^, M^ = m'p + nm'Q, n = l,2,...; Cq 
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and the corresponding Borel transform is 

2 



2 

^'^^^^ = 3M^Jo 



Substituting ( |lOlD into (|To|) yields 

2 ( n oo 
/o(M) = ^ \ -e-"p/^^' + ^ e-K+nm2)/M2 I 
I ^ n=l J 

^^^2m(^)\ ™S = 4,a, (103) 

vr \ vr / 

This should be compared to the standard result ^ with standard (small) condensates 



TT^ G2 ^ 448^K(teW. (104) 



3 M4 81 Af6 
In ( |104D as{M) is standard, i.e. obtained from by setting Mb = 0. 



2 

It is clear that (|103| ) contains in the Borel plane a set of poles at = M| = ±^2;^^, k = 
1,2,... and an essential singularity at M = 0. These features imply presence of renormalons and 
are connected to the factorial growth of coefficients A2n in (|100|). 



Now remarkably both Borel transforms lie inside the corridor of experimental errors, thus 
describing satisfactorily data with very different values of condensates (for details of comparison 
see [^2|). Thus situation is becoming even more unclear: not only one has two sets of condensates 
(and consequently two sets of sum rules) but in addition experimental data cannot give preference 
to one of them. 

While leading perturbative large-M asymptotics of loiM) and Iq{M) coincide, there is an 
important difference at small M: while M'^Iq{M) is defined for all M, M'^Iq^{M) is diverging 
for M — > due to higher condensates and higher powers of as{M). 



8 Conclusions and outlook 



The main emphasis of the present paper is the influence of conflnement on the behaviour of 
Green's functions in their dependence on momentum and the behaviour of Borel transforms. 
We stressed above everywhere the importance of large distances working in coordinate represen- 
tation, especially for light quarks in presence of conflnement. As a first and most clear example 
the Green's function of Dirac equation with linear scalar potential was considered and it was 
demonstrated that the Eucliden time expansion (equivalent of Borel transform for heavy-light 
systems) looks completely diff'erent from the nonrelativistic case, and from the template oscilla- 
tor Green's function. In this way it was shown that large distances may be important even for 
small Euclidean times and bring about new terms in the OPE in coordinate space. 



21 



As a second example we have treated the nonhnear equations for a quark in the heavy-hght 
system -nonlocal equivalent of the Dirac case, and found that again the result is different from 
what one would expect in standard OPE, but the terms of expansion turn out to be constant. 
Si (T) ~ const ■ a^l'^ . Translating this contribution into the form of the usual correlation function 
n((5^) of vector currents (like it is done in the reaction e^e~ — > hadrons) one would have the 
contribution An((5^) ~ , which is similar to the standard term -^P^, and is presumably 
one term in the subseries generating ra{qq). In this example large distances, explicitly accounting 
for in our analysis, do not produce new OPE terms but give some path to calculating chiral 
condensate through confinement characteristics (i.e. string tension a). 

In section 4, in contrast to that, another problem was elucidated: how linear confinement 
is built up out of higher condensates of OPE, and the answer is given by comparison of Eqs. 



and (51). Indeed, the infinite sum of derivatives of field correlators in ( |5l|) is equivalent 
to the linear confinement term in (^9[), and to extract it explicitly one needs to rearrange all 
derivatives. 

We have analyzed abelian confining models in section 6 and described different possible 
sources of nonstandard OPE terms, e.g. ■ 

Finally, the last problem considered in the paper concerns the derivation of OPE from 
the spectral representation of the meson Green function. When the spectrum and coefficients 
Cn (equivalent of quark decay constants /tt) are known, the OPE is calculated automatically 
and can be compared with that obtained "microscopically"- i.e. via Feynman diagrams in the 
external fields and equations of motion. 

In the d = 1 + 1 QCD this program was fully investigated in a series of papers (see e.g. [^] 
and refs. therein) and a mismatch between condensates obtained in those two ways was found. 



In the 3+1 QCD situation is similar and as shown in |32| and in the present paper, the 
mismatch of condensates in scales and order of magnitude also is evident. The situation is 
sharpened by the fact, that the QCD sum rules for e^e~ hadrons reproduce experimental 
data for both choices of condensates. 

We have not tried here to resolve this puzzle, and leave it for the future. There are two 
important topics in OPE we have not discussed. First, this is the partial summation of the OPE 
terms which can be done by introduction of nonlocal condensates in OPE, initiated and studied 
in |p7|, 38 1 . It would be interesting to find the link between our treatment of long-distance 



nonperturbative physics and the method of nonlocal condensates worked out in |37 



Second, the problem of perturbative-nonperturbative interference, which may produce new 
singular OPE terms, like 1/Q^, which was discussed in |4^, |9|, is touched in section 6 only 
briefly. This set of problems certainly deserves futher study. 
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A Appendix 



We discuss in this Appendix the properties of the kernel (18), which we used in the main text. 
The reader is referred to the Appendix 3 of the paper where 3d counterpart of (|l^) was 



analysed. Since the kernels of the form (18) play important role in the discussed formalism 
we present an independent detailed analysis here both for possible future applications and for 
reader's convenience. 

We are interested in the properties of the following function: 

1 1 , ^ ^ 2\ 

f{ff,p) = jdaaj d/3/3exp (- ^""^"/^^ J (A.l) 

\ 9 / 

where fj, p are d-dimensional vectors with angle between them. We denote absolute values of 
the arguments as r] = \ff\, p = \f)\ and assume in what follows without loss of generality that 
r] > p. The symmetry of formulas below with respect to the exchange p ^ r] (which is manifest 
in the definition ( A.l] )) is to be restored by replacements p min {/?, ry} and r] max {p, rj}. 



It is instructive to consider four different asymptotic regions: 

l).r],pr^Tg 2). ri^^Tg-p^Ty 3). r?, p » T,,; > 1 4). ??, p » T^; « 1 



In the region 1) one can expand ( [A.l| ) in Taylor series with respect to both arguments, 
subsequent integration is straightforward: 

^ 1 p^ + if , 2prjcose , p^ + r/^ , p^r?^ f 2 n , 

/(^>A) = 4-^^ + ^^ + ^i^ + ^(cos e + - 



2p7] cos 6 



(p2 + ^2)+0(r-6) (A.2) 



In derivation of expression (^) in the main text we have used in fact the leading term of this 
asymptotics (i.e. 1/4). 

In the regions 2), 3), 4) we will systematically omit exponentially small terms, i.e. terms 
proportional to exp(— r/^/T^) and also terms ~ exp(-pV^g ) in the regions 3) and 4). One 
easily obtains the following expression in the region 2) : 

Now we come to the regions 3) and 4). It is instructive to introduce the following variables: 

{fjp - pr]f Aifp'^ .^9 {ffp + pr])^ rfp^ ^9 
s = 7^ = ^„ sm - ; q= — = ^„ cos - (A. 4) 

r2 r| 2 ' ^ 4T| r2 2 ^ ^ 

and also ^ = a/s / = 2 tan | . 
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In the region 3) the upper hmit of the integration in ( |A.1| ) can be shifted to infinity up to 
exponentiahy smah corrections. The function f{fj,p) in the region 3) can be written therefore 

as ^ 



where (p{(,) is given by 



x/^ (4 + e 
8 



2^2 



J (iyexp(-y2) 



1 - Erf ( ^ 



(A.6) 

The function (/>(^) is a monotonically decreasing function of ^. When ^ is going to infinity, 0(^) 
is approaching the following asymptotics: 



1 "^9 1 

*(0 = 3 + i^|j+o(r 



(A.7) 



At the point = 2, which corresponds io 6 = tt/2 and hence orthogonal vectors r/ and p one 
finds (/>(2) = 1, in agreement with simple direct calculation from (|A.lD . 

We are now in the position to analyse the properties of /(??, in the region 4), where 
^ ~ ^ 1. Making the change of variables, one gets from ( [A.l[ ): 



fi'n,p) 



rp2 

a 

sm.0 



V 2 



PC P ^ 2 



-w 77 2 



y J dx + J dy J dx + J dy J dx— 



Ml 
2 



Ml 

2 



Ml 

2 



p ^ 2 



dy 



dx^\^- y"^" 



4s 



exp i^—x"^ — 



2 



(A. 



where w = yjq{l/ p — l/i]) = (rj — p) cos (0/2) /T^. One can rewrite ( |A.8| ) in the following form 



fif],p) 



9 ^ 



smf s 



V 



(A. 



where the ^-expansion of the functions g, /s, /4 can be performed systematically. It gives 



vr 



exp(— 



vr 



2^3 



C--Uiz) + —exp{-z')].^' + 0{e 



vr 



/2(z, t., 6 = Erf {w)^k{z) + [1 - exp(-u;2)]^^PL£V . ^ + 



(A.IO) 
(A.ll) 



where the function k(z) is defined as 



k{z) =Erf(z) 



2z 



exp(— z^) 



(A.12) 
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Extracting coefficient functions in front of higher powers in ^ is a matter of straightforward 
algebra. 



The expressions ( A.IOD , ( A.ll ) are exact at the given order in ^ up to omitted exponentially 
small terms. They can be simplified in different limiting cases. If w = (i.e. tj = p ), one has 
/2 = 0, while the first two terms in the r.h.s. of ( A.9| ) are equal. In the opposite limit w ^ oo 
the following relations hold true: 



lim 




lim /2 



(A.13) 



Notice that in all cases the first argument of the functions g, f2 need not be small: y/s/r] 
{2p/Tg) sm{e/2) and in the region 4) /? > Tg, but 6* < 1. 

In terms of the original variables rj, p, 9 the leading term in ( |A.9D can be represented as 



f{v,p) 



Z9_ 

64 sin^ I cos | rfp'^ 



2p sin ■ 



(1 + Erf(u')) +K 



2r/ sin ■ 



;i - Erf (li;)) 



f / \ 9 / 

(A.14) 

where w = (rj — p) cos{6 /2)/Tg and k{z) is defined in ( A.12| ), k{z) > if z > 0. This expression 
is valid in small 0-limit. 



Notice that / is non-singular if ^ limit (which is evident from ( |A.l| )): 

Tg^ 



4 + 4+Erfte^^^^ ^ 



jj2 p2 



(A.15) 



One needs some simple extrapolating representation of (A.l) for practical calculations. 
Notice that it is only asymptotic behaviour of /(??, /o) that matters, the particular form of 
Gaussian kernel was taken in ( [A.l| ) just as an example. A possible expression respecting all 
desired properties of / in the regions of large r/, p is as follows: 



f{ff,p) 



rp4 



4fj2p 



(A.16) 



where the function 1(9) has the following "focusing" property: being integrated with a regular 
function F{9), it acts like a smoothed 5-function (see p^): 



d9F{9)m ^ c,^^ F ('-^ 



+ C2J d9F{9) 



(A.17) 



where cq,ci,C2 are some constants of the order of unity. It is worth reminding that p is the 
length of the smaller vector in our notation, i.e. p = min{p, r/}. In particular, it is seen that 
in the limit of large p ^ Tg the small 9 asymptotics gives dominant contribution unless F{9) 
vanishes at the origin faster than 9"^. 
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